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Abstract 

The  onset  of  axisymmetric  (bubble)  and  non-axisymmetric  (spiral)  modes  of  breakdown  is  stud¬ 
ied  numerically  for  swirling  pipe  flows.  We  find  that  the  onset  of  axisymmetric  vortex  breakdown 
occurs  when  the  vortex  attains  local  criticality.  A  transient  simulation  of  the  evolution  of  vortex 
breakdown  revealed  that  downstream-running  waves  are  trapped  approximately  at  the  location 
of  flow  criticality.  These  trapped  waves  are  slowly  amplified  and  eventually  result  in  the  bubble 
breakdown  region  with  reversed  flow  and  enlarged  core  size.  Non-axisymmetric  disturbance  were 
found  to  decay  on  columnar  base  flows;  however,  for  base  flows  with  bubble  breakdown,  non- 
axisymmetric  disturbances  were  amplified.  These  amplified  disturbances  resulted  in  the  formation 
of  spiral  breakdown.  If  the  base  flow  inlet  swirl  was  only  slightly  larger  than  that  leading  to  ax¬ 
isymmetric  breakdown,  a  nearly  columnar  solution  was  obtained,  suggesting  that  weak  asymmetry 
may  help  to  stabilize  the  columnar  solution  branch.  The  numerical  simulations  also  revealed  that 
most  of  the  energy  in  spiral  breakdown  is  contained  in  the  first  few  non-axisymmetric  modes  of 
variation. 


1  Introduction 

Since  first  observed  by  Peckham  and  Atkins [1],  the  sudden  transition  of  a  concentrated  vortex- 
dominated  flow  to  a  less  coherent  spiraling  flow  referred  to  as  vortex  breakdown  has  attracted  consid¬ 
erable  attention  for  both  practical  and  theoretical  reasons.  The  appearance  of  vortex  breakdown  over 
swept  wings  at  high  angles  of  attack  imposes  additional  requirements  on  manuevering  for  these  vehi¬ 
cles  and  can  lead  to  significant  design  issues  associated  with  increased  fatigue  on  parts  of  an  aircraft 
in  its  wake.  The  use  of  the  recirculation  region  associated  with  vortex  breakdown  as  a  fluid  dynamic 
flame  holder  in  combustion  engines  provides  a  mechanism  to  increase  the  mixing  between  the  fuel  and 
the  air  without  needing  to  decelerate  the  main  flow  as  significantly  as  would  otherwise  be  necessary, 
thereby  increasing  the  efficiency  of  the  engine. 

Although  there  are  two  dominant  modes  for  vortex  breakdown,  as  many  as  seven  different  forms 
have  been  cataloged  by  Faler  and  Liebovich[2]  and  by  Sarpkaya[3].  The  “bubble”  mode  is  predomi¬ 
nantly  axisymmetric  in  nature  and  is  characterized  by  a  recirculation  region  several  times  larger  than 
the  vortex  core  prior  to  breakdown.  The  “spiral”  mode  tends  to  be  observed  at  lower  swirl  ratios  than 
the  bubble  mode  and  is  distinguished  by  a  sharp  bend  in  the  vortex  core,  followed  by  the  formation 
of  several  turns  of  a  helix  before  losing  its  coherency.  The  different  modes  are  sensitive  to  flow  condi¬ 
tions  and  have  often  been  observed  to  spontaneously  transition  from  one  mode  to  another.  Brucker[4] 
used  laser  sheet  visualization  to  more  closely  examine  the  transition  between  modes  and  identified  the 
development  of  a  rotating  structure  inside  the  recirculation  zone  of  a  bubble  mode  that  tilted  and 
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extended  to  create  the  helix  of  the  spiral  mode  as  the  swirl  ratio  was  decreased.  The  reverse  of  this 
procedure  occurred  when  the  spiral  mode  transitioned  to  the  bubble  mode. 

The  main  emphasis  of  this  paper  is  to  explore  the  connection  between  wave  propagation  and  the 
onset  of  bubble  and  spiral  modes  of  vortex  breakdown.  Squire[5]  is  generally  credited  with  the  first 
suggestion  that  vortex  breakdown  could  be  explained  with  concepts  of  wave  motion.  A  central  element 
of  Squire’s  work  (and  most  wave  theories  of  vortex  breakdown  since  Squire)  is  the  division  of  vortical 
flows  into  supercritical  and  subcritical  states.  Supercritical  vortices  support  only  downstream-running 
waves  while  subcritical  vortices  support  both  downstream-  and  upstream-running  waves.  Experimental 
data  shows  that  vortices  without  breakdown  or  upstream  of  breakdown  are  supercritical[6,  7,  8].  Squire 
proposed  that  vortex  breakdown  becomes  possible  when  a  vortex  first  achieves  criticality.  As  discussed 
in  numerous  reviews[8,  9,  10,  11,  12],  Squire  suggested  that  disturbances  that  are  present  downstream 
could  be  transmitted  upstream  when  a  vortex  achieves  critical  conditions,  and  these  upstream-running 
disturbances  would  then  be  responsible  for  initiating  vortex  breakdown.  Thus,  according  to  Squire’s 
hypothesis,  vortical  flows  at  the  onset  of  breakdown  should  attain  critical  conditions  at  some  location 
along  the  vortex  core.  Benjamin[13,  14]  advocated  a  similar  argument  based  on  the  existence  of 
conjugate  flow  states  at  the  breakdown  location. 

Recent  theoretical  and  computational  work  has  resulted  in  significant  progress  in  the  understanding 
of  axisymmetric  breakdown.  Randall  and  Leibovich[15]  used  a  weakly  nonlinear  analysis  to  show  that 
the  steady,  Euler  equations  contain  a  solution  branch  which  bifurcates  from  the  columnar  vortex 
solution  at  the  critical  swirl.  This  branch  of  solutions  is  characterized  by  a  standing  solitary  wave 
which  evolves  from  a  columnar,  swirling  base  flow  in  an  infinitely  long,  straight  pipe.  Leibovich  and 
Kribus[16]  have  extended  this  to  larger  amplitudes  and  show  that  a  stagnation  point  can  develop  in 
the  solitary  wave. 

Beran  and  Culick[17]  have  solved  the  steady,  axisymmetric  Navier-Stokes  equations  in  pipes  of 
varying-area  and  revealed  a  more  complicated  bifurcation  picture.  Specifically,  for  sufficiently  high 
Reynolds  number,  they  found  the  existence  of  two  limit  points.  As  shown  in  Fig.  1,  limit  point  A 
marks  the  end  of  the  near  columnar  solution  branch  and  is  the  highest  value  of  swirl  for  which  a 
columnar  steady  solution  was  found.  Above  this  swirl,  the  solution  jumps  to  the  lower  branch  which 
contains  bubble  breakdown  with  large  regions  of  reversed  flow.  Limit  point  B  marks  the  end  of  the 
breakdown  solution  branch.  For  lower  swirls,  the  solutions  must  return  to  the  columnar  state.  In 
between  the  two  limit  points,  the  solutions  contain  a  localized  standing  wave  which  gradually  merges 
with  the  other  branches.  More  recent  analysis  has  shown  this  localized  standing  wave  branch  to  be 
unstable[18,  19]. 

Although  Beran  and  Culick  did  not  possess  definitive  proof,  they  suggested  that  limit  point  A 
corresponds  to  the  occurrence  of  critical  conditions  in  the  flow.  Specifically,  they  solved  Hall’s[10,  20] 
quasi-cylindrical  equations  and  showed  that  numerical  integration  of  these  equations  fail  at  approxi¬ 
mately  the  same  inlet  swirl  as  limit  point  A  in  the  full  Navier-Stokes  simulations.  Since  Hall[10]  has 
shown  that  the  quasi-cylindrical  equations  fail  at  the  critical  state,  Beran  and  Culick  hypothesized 
that  limit  point  A  corresponds  to  the  occurrence  of  criticality  in  the  flow,  and,  thus,  is  equivalent  to 
the  bifurcation  described  by  Randall  and  Leibovich.  One  contribution  of  this  paper  is  to  show  that 
this  limit  point  does  indeed  coincide  with  the  occurrence  of  local  flow  criticality. 

In  a  recent  series  of  papers  by  Wang  and  Rusak[21,  22,  23,  24],  significant  advances  have  been  made 
to  unify,  in  a  solid  theoretical  foundation,  the  various  studies  of  axisymmetric  vortex  breakdown.  Their 
basic  approach  is  to  study  the  stability  of  columnar  solutions  of  the  Euler  equations  in  a  finite  length, 
constant  area  pipe  and  connect  the  stability  results  to  the  time  asymptotic  solutions  supported  in 
the  pipe.  In  their  theoretical  analysis,  they  find  a  bifurcation  diagram  similar  to  that  of  Beran  and 
Culick’s  numerical  study.  According  to  their  theory,  limit  point  A  occurs  at  the  equivalent  of  Squire’s 
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A 


Figure  1:  Minimum  axial  velocity,  for  steady  solutions  of  the  Navier-Stokes  equations  versus  inlet 
swirl,  $7,  as  observed  by  Beran  and  Culick[17]. 


or  Benjamin’s  critical  state  for  a  finite  length  pipe.  Above  the  critical  swirl,  the  columnar  solutions 
are  unstable  and  the  breakdown  branch  must  be  observed  as  the  time-asymptotic  state.  Limit  point 
B  is  equivalent  to  Escudier  and  Keller’s  two-stage  transition  model[25]  of  vortex  breakdown  which 
features  a  hollow  breakdown  region  with  stagnant  flow.  In  numerical  simulations  of  the  unsteady, 
Euler  equations [26],  they  confirm  their  analysis  and  show  that  finite  amplitude  disturbances  can  result 
in  a  columnar  solution  developing  into  a  breakdown  solution. 

The  remainder  of  the  paper  proceeds  as  follows.  First,  we  concentrate  on  axisymmetric  pipe 
flows  and  use  numerical  simulations  to  investigate  the  onset  of  bubble  breakdown.  In  particular,  we 
show  that  the  onset  of  bubble  breakdown  coincides  with  local  flow  criticality.  Then,  we  consider 
non-axisymmetric  disturbances  on  axisymmetric  base  flows  and  use  transient  simulations  to  study 
the  onset  of  spiral  breakdown.  The  simulations  show  that  the  non-axisymmetric  disturbances  do  not 
grow  until  bubble  breakdown  is  present  in  the  axisymmetric  base  flow.  In  addition,  we  show  that 
the  majority  of  the  energy  content  of  spiral  breakdown  is  contained  in  the  first  few  non-axisymmetric 
modes. 

2  Axisymmetric  Flow 

We  begin  by  analyzing  the  criticality  of  ten  different  pipe  flow  trials  which  have  been  simulated  using 
an  axisymmetric  Navier-Stokes  solver.  The  algorithm  used  to  calculate  the  incompressible,  swirling 
flow  is  a  finite  volume  discretization  described  in  detail  in  Darmofal[27].  The  solver  has  been  validated 
with  numerous  comparisons  to  both  experiments  and  computations  by  other  investigators [2 7].  Of 
specific  interest  to  this  work,  the  solver  has  been  used  to  show  that  axisymmetric  simulations  provide 
good  quantitative  agreement  with  experimental  results  for  bubble  breakdown  in  swirling  pipe  flows[28]. 

A  single  trial  for  a  given  pipe  flow  consists  of  starting  with  a  zero  inlet  swirl  steady  flow  and  quasi- 
statically  raising  the  inlet  swirl  until  breakdown  occurs  at  a  swirl  of  Qf,.  The  criticality  is  determined 
by  solving  an  appropriate  eigenvalue  problem  as  described  by  Benjamin [13]  and  in  a  more  approximate 
manner  using  a  local  mean  swirl  number  similar  to  Spall[29],  Robinson,  et  al.  [30],  and  Delery[12],  The 
eigenvalue  criticality  calculation  assumes  the  flow  can  be  analyzed  using  locally  periodic  eigenmodes 
and  thus  requires  small  axial  gradients.  Since  the  eigenvalue  analysis  is  conducted  on  only  smooth, 
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R0 
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25 

30 

2.0 

1.8 

2.0 

1.8 

B 

5 

10 

25 

30 

2.0 

1.8 

2.4 

1.8 

Pipe 

Zi 

Zd 

Zt 

•^max 

Ri 

Rt 

Rc 

C 

3.4 

70 

75 

80 

5.0 

6.67 

5.0 

D 

3.4 

37 

42 

47 

5.0 

5.84 

5.0 

Table  1:  Constants  for  pipe  geometries 


nearly  columnar  flows,  these  conditions  can  generally  be  satisfied.  Our  results  show  conclusively  that 
the  limiting  swirl  at  which  breakdown  occurs  coincides  with  the  occurrence  of  local  flow  criticality  as 
suggested  by  both  Beran  and  Culick  and  by  Wang  and  Rusak. 

The  inlet  boundary  conditions  model  the  incoming  flow  as  a  g-vortex[8,  31].  Specifically, 

u(r,  z  =  0)  =  0, 

F(r,z  =  0)  =  n  [1  -  exp(-r2)]  ,  (1) 

w(r,z  =  0)  =  1  +  Aw  exp(—  r2). 

Velocities  are  non-dimensionalized  by  Woo,  the  free  stream  axial  velocity.  The  circulation  function  is 
defined  as  T  =  rv  and  the  value  of  q  as  defined  by  Batchelor[31]  is  q  =  £1/ Aw.  The  vortex  core  size, 
<5,  is  used  as  a  length  scale.  Therefore,  r  =  1  is  the  core  edge,  even  though  the  maximum  swirl  occurs 
at  r  =  1.21.  The  pipe  wall  is  modeled  as  an  inviscid,  slip  boundary  in  order  to  avoid  any  ambiguities 
associated  with  wall  boundary  layer  behavior.  Finally,  at  the  outlet,  the  flow  is  assumed  to  have 
negligible  axial  gradients  and  all  axial  derivatives  are  approximated  by  one-sided  differences.  For  all  of 
the  axisymmetric  results  in  this  paper,  the  Reynolds  number,  Re  =  W^S/u,  is  set  to  1000.  We  have 
found  no  substantial  variations  in  flow  characteristics  for  Reynolds  numbers  above  about  500[27]. 

Two  different  pipe  geometries  have  been  studied  in  this  paper.  The  first  possesses  a  converging- 
diverging  inlet  section,  a  constant  area  test  section,  and  converging  outlet  section;  it  is  a  slight  variant 
of  the  geometry  considered  by  Beran  and  Culick[17].  The  inlet  section  serves  to  isolate  the  inlet 
boundary  from  the  breakdown  bubble  which  occurs  in  the  diverging  section  of  the  pipe.  The  specific 
pipe  radius,  R(z),  chosen  for  this  investigation  is 

Ri  +  {Rt  -  Ri)g{z,zt ), 

Rt  +  {Ro  -  Rt)g{z  -  zt,  z0  -  zt), 

Ro, 

Ro  T  (Re  Ro)g(z  zc,  zmax  zc), 

g(zj)  =  ^{1  -  cos[7r(*/0]}  • 

Specific  constants  for  the  two  pipes  (A  and  B)  with  this  basic  form  are  given  in  Table  1.  Note  that 
Pipe  B  has  a  larger  area  increase  than  Pipe  A. 

The  second  geometry  was  used  in  the  experiments  of  Sarpkaya[3,  32,  33]  and  Leibovich  et  al.[ 2,  6, 
34],  This  pipe  consists  of  a  constant  area  inlet  of  radius  Rj  and  length  z,L.  a  divergence  section  from 
Zi  <  z  <  Zd  with  a  linear  pipe  radius  variation,  a  constant  area  test  section  from  z([  <  z  <  zt  with  radius 


<  z  <  zt 

<  z  <  zQ 

<  Z  <  Zr 


(2) 
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Pipe 

Fixed  q  value 

fib 

5 

A 

q  =  1.8 

3.91 

6 

B 

q  =  1.6 

2.92 

7 

B 

q  =  1.8 

2.41 

8 

B 

q  =  2.0 

2.10 

9 

C 

q  =  1.8 

1.66 

10 

D 

q  =  1.8 

2.15 

Trial 

Pipe 

Fixed  Aw  value 

1 

A 

Aw  =  0.0 

1.51 

2 

A 

Aw  =  0.4 

1.97 

3 

B 

Aw  =  0.0 

0.91 

4 

B 

Aw  =  0.4 

1.38 

Table  2:  Trial  definitions  and  values 


Rt,  and  a  converging  outlet  section  from  zt  <  z  <  zm ax.  The  specific  form  of  the  converging  outlet 
section  was  not  given  in  the  references  so  we  employ  the  same  sinusoidal  variation  used  in  the  outlet 
sections  of  Pipes  A  and  B.  We  consider  two  pipes  (C  and  D)  with  specific  constants  given  in  Table  1. 
Both  pipes  have  the  same  divergence  angle  as  that  used  in  the  previous  experiments.  Pipe  C  has  the 
full  length  of  the  divergent  section  used  in  the  experiments  while  Pipe  D  has  approximately  half  the 
length.  As  we  show  below,  while  the  quantitative  details  are  different  for  these  pipes,  the  qualitative 
results  are  unchanged.  Also,  both  pipes  feature  a  constant  area  test  section  which  is  significantly 
shorter  than  that  in  the  experiment;  again,  this  change  should  only  affect  the  quantitative  details. 

The  inlet  radius  of  both  Pipes  C  and  D  is  5.(W  which,  when  using  the  (/-vortex  for  the  inlet  condition, 
implies  the  maximum  swirl  occurs  at  0.242it!j.  This  was  chosen  to  provide  a  realistic  comparison  with 
the  previous  experiments.  For  example,  Leibovich  and  Faler[2]  present  velocity  profiles  with  the 
location  of  maximum  swirl  ranging  from  approximately  1/4  to  1/3  of  the  pipe  radius. 

In  all  simulations,  the  axial  grid  spacing  was  A z  =  0.2.  For  Pipes  A  and  B,  the  radial  distribution 
of  points,  rj,  was  constant  at  a  given  axial  location, 

rj  =  R(z)  J  *  , 

Jmax  -L 

with  jmax  =  31.  The  radial  distribution  of  points  for  Pipes  C  and  D  was  designed  to  decrease  resolution 
outside  the  region  where  the  core  was  expected  and  satisfied, 


for  1  <  j  <  jc ,  and 


rj  =  fcR(z)~ 


3  ~  1 


jc  -  1 


Tj  =  R(z) 


fc  +  (1  “  fc ) 


3  ~  Jc 

Jmax  —  jc 


for  jr  <  j  <  jmax-  Specifically,  we  used  fc  =  0.4,  jc  =  31,  jmax.  —  51,  and  c  =  1.27.  The  cluster  param¬ 
eter,  c,  was  chosen  such  that  the  first  cell  outside  of  the  equally-spaced  core  region  was  approximately 
the  same  size  as  the  cells  in  the  core  region. 

As  described  above,  ten  trials  were  run  varying  the  inlet  swirl,  S7,  in  a  quasi-static  manner.  In 
these  trials,  either  the  value  of  Aw  or  q  is  held  constant  at  the  inlet.  The  first  solution  of  a  trial  is 
found  by  starting  from  an  undisturbed  flow  with  zero  swirl.  After  finding  the  steady  state  solution, 
Q  is  raised  and  the  solution  is  again  converged.  The  trial  is  complete  when  breakdown  first  occurs  at 
the  onset  swirl,  Vtb.  The  value  of  Qb  is  determined  to  the  nearest  hundredth.  The  specifics  of  the  ten 
trials  described  in  this  paper  and  the  resulting  onset  swirls  are  presented  in  Table  2. 
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(a)  Q  =  0.50 


w  G  (0.9, 1.3),  Aw  =  0.025 


Streamsurfaces  w  G  (0.9, 1.4),  Aw  =  0.025 

(b)  n  =  i.oo 


Streamsurfaces  w  G  (0.5, 1.6),  Aw  =  0.05 

(c)  n  =  1.50 


Streamsurfaces  w  G  (—0.2, 1.6),  Aw  =  0.10 

(d)  n  =  1.51 


Figure  2:  Streamsurfaces  and  axial  velocity  contours  for  Trial  1  axisymmetric  steady  solutions.  Re  = 
1000. 


Sample  streamsurfaces  and  axial  velocity  contours  from  four  steady  solutions  of  Trial  1  are  shown  in 
Figure  2.  For  this  trial,  f \  =  1.51.  The  streamsurfaces  and  the  axial  velocity  contours  for  the  solutions 
with  12  <  SI;,  show  no  indication  of  the  impending  breakdown  at  12  =  12;,.  Although  the  minimum  axial 
velocity  is  decreasing  with  increasing  swirl,  the  flow  is  not  near  stagnation  for  12  =  1.50.  An  important 
point  to  notice  is  that  even  the  flows  which  are  on  the  verge  of  breakdown,  such  as  12  =  1.50,  are  still 
smoothly  varying.  This  provides  some  assurance  that  the  criticality  analysis  in  the  following  section 
is  not  flawed  by  the  presence  of  large  axial  variations.  The  smoothness  of  the  flows  is  especially  true 
in  the  constant  area  section  of  the  pipe  which  will  be  the  first  location  to  achieve  critical  conditions. 
Finally,  as  discussed  by  Beran  and  Culick[17],  the  drastic  change  in  solution  behavior  is  a  result  of 
the  solution  reaching  a  limit  point  in  the  solution  space  and  jumping  to  the  steady  solutions  with 
breakdown.  In  the  next  section,  we  will  show  that  the  limit  point  and  corresponding  jump  is  a  result 
of  the  vortex  becoming  locally  critical.  Although  we  do  not  present  any  other  streamsurface  and 
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velocity  plots,  this  sudden  onset  of  vortex  breakdown  was  characteristic  of  all  ten  trials. 


2.1  Axisymmetric  criticality  analysis 

For  axisymmetric  flows,  it  is  possible  to  determine  the  criticality  of  a  flow  using  a  standing  wave 
analysis[13].  To  analyze  standing  waves,  the  vortex  flow  is  divided  into  a  steady,  mean  flow  and  a  steady 
infinitesimal  perturbation.  The  steady,  mean  flow  velocity  distribution  is  given  by  [0,  V (r,  z),  W (r,  z)\ 
and  is  assumed  to  have  axial  variations  whose  length  scale,  L ,  is  large  compared  to  the  perturbation 
wavelength  such  that  |y| L  »  1.  The  streamfunction  perturbation  is  defined  as  ip  =  f(r)exp(,yz) 
where  7  is  the  axial  wavenumber.  As  given  by  Hall[10],  the  governing  equation  for  stationary  pertur¬ 
bations  is: 


r±  n  c£\ 

dr  \r  dr  J 


+ 


r  d  /ldffA 
W  dr  \r  dr  ) 


+ 


1  d(r2)~ 

r3W2  dr 


f  =  0, 


(3) 


where  T  =  rV(r,  z)  is  the  mean  flow  circulation.  For  flow  in  a  pipe,  the  perturbation  streamfunction  is 
zero  at  r  =  0  and  at  the  pipe  wall,  r  =  R(z).  Hall[10]  has  shown  that  Eq.  3  is  valid  for  quasi-cylindrical 
flows  as  well  as  columnar.  An  infinite,  ordered  set  of  eigenvalues  exists  for  this  Sturm-Liouville  system 
such  that  7q  <  y2  <  yf  <  ...[13].  Thus,  the  flow  is  supercritical  when  Jq  >  0  since  then  all  eigenvalues 
of  Eq.  3  will  be  positive  and  only  exponential  solutions  are  possible.  When  Jq  <  0,  some  perturbations 
are  sinusoidally  varying  waves  and  the  flow  is  subcritical.  In  the  supercritical  case,  I70 1  gives  the  slowest 
decay  rate  of  disturbances  in  the  upstream  direction  and  is  a  measure  for  how  much  a  vortex  will  be 
influenced  by  downstream  effects.  For  subcritical  flows,  1 70 1  is  the  highest  wavenumber  possible  for 
stationary  waves  on  the  mean  flow. 

In  the  following,  the  criticality  of  the  132  steady  solutions  from  the  10  trials  will  be  determined  by 
calculating  the  local  eigenvalues,  y2,  at  each  axial  grid  station  for  Eq.  3.  The  local  mean  flow  is  given 
by  the  local  velocity  distributions  from  the  steady  results  of  the  axisymmetric,  Navier-Stokes  solver. 
The  eigenvalues  are  calculated  by  discretizing  Eq.  3  with  a  second  order  accurate  finite  difference 
scheme  and  using  an  EISPACK  eigenvalue  solver  on  the  resulting  matrix.  For  display  purposes,  only 
the  value  of  Jq  will  be  plotted  since  this  eigenvalue  completely  determines  the  local  criticality. 

We  have  also  examined  some  approximate  criteria  which  are  often  used  to  determine  whether  or 
not  breakdown  is  impending.  Specifically,  an  increase  of  swirl  ratio  generally  results  in  a  decrease  of 
vortex  criticality  from  supercritical  towards  subcritical.  Numerous  authors[12,  29,  30]  have  investigated 
different  techniques  for  evaluating  the  swirl  ratio.  For  some  unique  swirl  ratio,  the  vortex  is  assumed 
to  become  critical  and  vortex  breakdown  occurs.  The  previous  studies  have  found  that  breakdown 
generally  occurs  for  swirl  numbers  above  approximately  1.4.  While  our  technique  for  determining  a 
local  swirl  ratio  is  again  different  from  these  references,  the  basic  ideas  are  the  same.  We  have  tested 
the  other  suggestions  for  evaluating  the  local  swirl  and  found  the  same  basic  trends  as  those  we  report 
using  our  method  although  our  method  seemed  to  provide  the  smallest  range  of  maximum  swirl  ratios 
at  breakdown  onset.  First,  we  define  the  local  swirl  number  as  S  =  Tc/Wcrc  where  rc,  Wc,  and  rc  are 
measures  of  the  local  circulation,  core  axial  velocity,  and  core  radius.  We  use  the  term  swirl  number 
instead  of  swirl  ratio  to  avoid  confusion  with  the  inlet  swirl  ratio,  ft.  The  core  radius  at  a  particular 
axial  station  is  determined  by  tracing  the  streamsurface  which  begins  at  the  inlet  location  of  maximum 
swirl  velocity.  For  a  given  z  location,  the  core  radius,  rc,  is  the  radial  position  of  this  streamsurface  at 
2.  The  axial  velocity,  Wc,  is  the  constant  axial  velocity  required  to  conserve  mass  in  the  core.  Finally, 
the  local  circulation  is  assumed  to  be  equal  to  the  inlet  value,  rc  =  ft. 

Figure  3  shows  the  axial  distribution  of  Jq,  the  centerline  axial  velocity,  w(0,z),  and  the  swirl 
number,  S  for  the  Trial  1  flows  with  12  =  1.48,  1.49,  1.50,  and  1.51.  Although  we  have  calculated 
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(a)  Minimum  eigenvalue,  7 q(z) 
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(c)  Swirl  number,  S(z) 

Figure  3:  7 q{z),  wq(0 ,z)  and  S(z)  distributions  for  Trial  1,  fl  =  1.48,  1.49,  1.50,  and  1.51,  Aw  =  0.0. 


the  eigenvalues  for  the  fl  =  fl/,  =  1.51,  they  should  not  be  interpreted  too  strongly  because  this  flow 
contains  breakdown  and  the  large  axial  gradients  surely  invalidate  the  eigenvalue  analysis  within  the 
vicinity  of  the  breakdown  region.  As  remarked  in  the  previous  section,  the  flows  with  fl  <  fl/,  are 
still  smoothly  varying  and  the  eigenvalue  analysis  should  be  reasonably  accurate.  The  results  show 
that  as  the  onset  inlet  swirl  ratio  is  approached,  the  vortex  is  becoming  locally  critical  at  z  ~  13. 
For  fl  =  1.50,  7q  is  nearly  zero.  By  contrast,  the  minimum  centerline  axial  velocity  is  approximately 
0.5  and  far  from  stagnation  for  fl  =  1.50.  This  highlights  that  vortex  breakdown  is  not  simply  a 
stagnation  process  but  rather  that  stagnation  is  a  consequence  of  vortex  breakdown. 

Next,  we  plot  the  variation  of  the  minimum  eigenvalue,  minimum  centerline  axial  velocity,  and 
maximum  swirl  number  versus  inlet  swirl  in  Figure  4  for  all  10  trials.  The  minimum  eigenvalue,  7om;nl 
and  minimum  centerline  axial  velocity,  ,  are  defined  as  the  minimum  value  of  these  variables  over 
all  axial  locations.  The  maximum  swirl  number,  5max,  is  similarly  defined.  Figure  4  is  a  compilation 
of  all  of  the  data  from  From  this  figure,  it  is  clear  that  the  criticality  as  determined  by  7om;n  is  a  good 
indicator  of  impending  vortex  breakdown.  The  minimum  eigenvalues  of  all  of  the  columnar  solutions 
with  fl  =  flfr  are  nearly  zero.  By  contrast,  the  swirl  number,  Smax,  is  not  as  precise;  breakdown 
occurs  for  values  of  5max  from  approximately  1.2  to  1.8.  While  the  swirl  number  may  not  be  useful  as 
a  precise  measure  for  the  onset  of  vortex  breakdown,  it  does  appear  to  give  the  correct  general  trends. 
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(a)  Minimum  eigenvalue,  7gmin  (b)  Minimum  axial  velocity, 


(c)  Maximum  swirl  number,  Smax 


Figure  4:  To  min,  Tnin  •  and  Smax  versus  54  for  all  steady  pipe  flow  solutions,  x,  columnar  solutions;  □, 
columnar  solutions  for  which  breakdown  occurs  if  54  is  raised  l/100th  (i.e.  54^  =  54  +  0.01). 


Specifically,  Smax  is  a  monotonically  increasing  function  as  breakdown  approaches;  thus,  any  effects 
which  raise  the  local  swirl  number  are  likely  to  drive  a  flow  closer  to  breakdown;  similarly,  effects 
which  lower  the  local  swirl  number  are  likely  to  stabilize  the  flow  away  from  vortex  breakdown. 

2.2  Transient  formation  of  axisymmetric  vortex  breakdown 

In  this  section,  we  report  the  results  of  a  transient  simulation  of  the  formation  of  vortex  breakdown. 
The  initial  flow  is  the  steady  result  from  Trial  1  for  54  =  1.50.  The  inlet  swirl  ratio  is  then  raised 
at  t  =  0  to  54  =  1.52,  initiating  a  disturbance  which  propagates  throughout  the  domain.  Since  this 
swirl  ratio  is  above  that  for  which  breakdown  occurs  (54 j  =  1.51),  the  final  solution  should  contain 
breakdown  and  the  transient  will  contain  the  breakdown  evolution.  Note,  the  solutions  in  this  section 
were  calculated  with  a  timestep  of  approximately  0.041  S/Wc 

Figure  5  shows  the  streamsurfaces  during  the  bubble  evolution.  At  t  =  41,  only  a  slight  per- 
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(a)  t  =  41  W^/5 
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c)  t  =  223  Woo  /S 

(h)  t  =  309  Woo/ 5 

(d)  t  =  254  Woo  5 


(i)  t  =  355  Woo  15 


(e)  t  =  286  Woo/5  (j)  Final  steady  result 

Figure  5:  Transient  bubble  formation  after  change  in  inlet  swirl  from  =  1.50  to  1.52  at  t  =  0+.  Pipe 
A,  Aw  =  0.0,  Re  =  1000. 
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turbation  is  visible  at  z  ~  12.4  By  t  =  165,  the  streamfunction  perturbation  has  increased  further, 
and,  at  t  =  223,  the  T  =  0  streamsurface  has  lifted  off  the  axis  to  form  a  small  region  of  reversed 
flow.  At  t  =  254,  the  bubble  has  become  quite  large  and  is  moving  upstream.  The  bubble  formation 
from  t  =  0  to  t  =  254  is  quite  symmetric  and  strongly  resembles  the  solitary  wave  solutions  found  by 
Leibovich[35],  suggesting  that  the  initial  development  of  breakdown  might  be  approximated  by  weakly 
nonlinear  theory.  At  t  =  286,  an  inner  bubble  is  formed  and  the  breakdown  bubble  loses  its  axial 
symmetry  while  continuing  its  upstream  motion.  As  the  breakdown  evolution  completes,  the  bubble 
takes  on  an  open-ended  shape,  and,  for  t  =  355,  a  second  bubble  forms  at  z  ps  17.6.  Finally,  the 
steady  result,  determined  when  the  residual  reached  machine  zero  at  approximately  t  =  1300,  shows 
the  first  bubble  has  increased  quite  significantly  in  size  and  the  second  bubble  has  moved  slightly 
upstream.  Evidently,  the  unsteady  bubble  may  be  smaller  than  the  final  steady  bubble  structure.  The 
calculation  was  continued  for  several  thousand  iterations  after  achieving  steady  state  with  no  changes 
in  the  solution  as  shown  in  Fig.  5(1). 

Escudier[9]  describes  an  experimental  breakdown  evolution  process  and  shows  a  time  sequence  of 
dye  visualization  photographs  which  are  remarkably  similar  to  the  previous  simulation  results.  In  the 
experimental  results,  a  slight  swelling  in  the  vortex  core  region  initially  appears  just  after  a  change 
in  flow  conditions.  This  swelling  gradually  develops  into  an  initially  axisymmetric  ring-like  structure. 
After  some  transient  asymmetries,  the  structure  stabilizes  into  a  bubble  breakdown  structure  slightly 
upstream  of  the  initial  swelling  location  with  a  second  stagnation  behind  the  primary  bubble. 

A  useful  technique  to  visualize  the  wave  motion  underlying  the  breakdown  evolution  is  to  track 
radial  perturbations  to  the  location  of  a  particular  streamfunction.  For  this  purpose,  we  track  the 
streamsurface,  tpma,Xi  which  begins  from  the  inlet  at  the  location  of  maximum  swirl.  Perturbations  to 
the  streamfunction  position  are  defined  implicitly  as, 

f{z,t)  =r{'ipmax,z,t)  -  r{ipmax,z,0). 

Thus,  regions  of  positive  r  represent  area  for  which  the  vortex  core  is  growing.  Figure  6  shows  the 
evolution  of  f  after  the  increase  in  the  inlet  swirl.  In  part  (a)  of  the  figure,  the  initial  stages  of  the 
development  are  shown.  A  small  amplitude  area- varying  wave  is  emitted  from  the  inlet  and  begins 
to  propagate  downstream.  A  portion  of  the  disturbance  propagates  across  the  domain  to  the  outlet; 
however,  some  of  the  disturbance  is  halted  very  near  the  location  of  flow  criticality.  As  shown  in  parts 
(b)  and  (c)  of  the  figure,  this  trapped  wave  continually  amplifies  and  slowly  moves  upstream  to  its 
final  position.  Thus,  one  might  interpret  this  breakdown  formation  as  resulting  from  trapping  and 
amplification  of  downstream-running  disturbances  at  the  critical  location. 

3  Asymmetric  Disturbances 

The  assumption  of  axisymmetry  provides  a  convenient  simplification  for  theoretical  development,  but 
since  physical  flows  are  not  bound  by  this  constraint,  theories  must  also  be  examined  in  the  absence 
of  this  simplification.  This  section  examines  the  behavior  of  a  weak  asymmetric  disturbance  on  a 
steady  axisymmetric  base  flow  characterized  by  a  given  inlet  swirl.  We  will  be  demonstrate  that, 
until  the  base  flow  reaches  a  particular  swirl  ratio,  the  disturbance  is  passively  convected  and  decays 
with  minimal  influence  on  the  base  flow.  However,  once  the  base  flow  achieves  a  critical  state  that 
would  exhibit  axisymmetric  vortex  breakdown,  significant  disturbance  growth  is  predicted  and  large 
perturbations  appear  on  the  base  flow.  The  basic  mechanism  for  this  growth  is  identified  through 
consideration  of  the  magnitude  of  the  velocity  components  for  each  Fourier  mode. 

A  psuedo-spectral  method  is  used  to  simulate  the  flow  by  solving  for  the  azimuthal  Fourier  modes 
of  the  velocity.  This  approach  is  more  efficient  for  computing  the  lower  azimuthal  modes  that  are  dom- 
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(a)  t  =  0  to  41  S/Woc  (b)  t  =  41  to  205  S/Woo 


(c)  t  =  205  to  533  S /Woo 


Figure  6:  Transient  evolution  of  the  radial  perturbation,  f,  of  the  maximum  swirl  streamline,  i/;max 
after  raising  inlet  swirl  from  52  =  1.50  to  1.52. 


inant  in  experiments  and  parallels  theoretical  techniques  for  investigating  the  flow.  Furthermore,  by 
using  the  Fourier-decomposed  model,  an  axisymmetric  solution  can  be  maintained  until  asymmetries 
are  deliberately  introduced.  The  Fourier-decomposition  of  the  velocity  field  is  specified  by 

Ne 

u  =  U  (r,z)  +u  o(r,z,t)  +  ^  (u  m(r,z,t)e  im0  +  um(r,  z,t)eim6j  ,  (4) 

m= 1 

where  U(r,  z)  is  a  steady  axisymmetric  base  solution  and  an  overbar  (”)  denotes  complex  conjugation. 
Thus  uTO  =  umer  +  vmeg  +  wmez  is  a  complex  function  of  r  z,  and  t  representing  the  mth  azimuthal 
mode  of  the  perturbation  velocity.  The  continuity  and  unsteady  incompressible  Navier-Stokes  equa¬ 
tions  are  written  in  cylindrical  coordinates  and  Fourier-decomposed  in  the  azimuthal  direction  to 
obtain  a  sequence  of  coupled  partial  differential  equations.  The  system  is  discretized  through  the  use 
of  Chebyshev  collocation  in  the  radial  direction  and  finite-differencing  in  the  axial  direction.  The 
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Figure  7:  Radial  distribution  of  axial  velocity  perturbation  mode. 


solution  is  obtained  using  an  operator  split  approach  that  successively  integrates  the  flow  in  time  by 
performing  a  non-linear  advection  step,  a  pressure  correction,  and  a  viscous  correction[36,  37]. 

As  in  the  axisymmetric  studies  from  the  previous  section,  the  flow  is  simulated  in  a  converging- 
diverging  nozzle  with  a  slip  wall  boundary  condition  on  the  nozzle  to  prevent  the  interaction  of  the 
boundary  layer  with  the  breakdown.  The  inlet  axisymmetric  velocity  distribution  is  specified  using 
Eq.  1.  Axis  boundary  conditions  consistent  with  the  Fourier-decomposition  are  used  and  the  outflow 
is  again  parabolized  using  one-sided  differences.  A  family  of  axisymmetric  solutions  are  obtained  using 
a  modified  Pipe  A  geometry  that  has  had  both  the  inlet  section  and  the  test  section  extended  by  5 
core  radii  and  the  final  converging  section  removed.  These  changes  were  made  to  allow  for  the  greater 
spatial  extent  of  the  anticipated  solution.  This  family  of  solutions  is  similar  to  that  obtained  in  Trial  1, 
but  was  computed  at  a  lower  Reynolds  number  of  500.  Despite  the  significantly  different  numerical 
scheme,  the  axisymmetric  solutions  with  the  pseudo-spectral  code  showed  very  good  comparison  to 
results  from  the  previous  axisymmetric  code[37]. 

Beginning  with  this  family  of  axisymmetric  base  flows,  a  weak  disturbance  of  the  form 

w\(r,z  =  0,  t)  =  A(t)r2e~12r  (5) 

is  introduced  to  the  m  =  1  mode  of  the  axial  velocity,  where  A  is  a  measure  of  the  amplitude  that  is 
linearly  ramped  in  time  to  1  percent  of  the  inlet  axial  velocity.  This  perturbation,  depicted  in  Fig.  7, 
is  chosen  to  satisfy  the  boundary  condition  at  the  axis  and  to  decay  rapidly  outside  the  vortex  core. 
Amplitude  values  ranging  from  0.05  to  5  percent  of  the  inlet  axial  velocity  have  been  simulated  with 
similar  qualitative  results  as  those  described  below.  As  a  consequence  of  the  physical  requirement 
for  single-valued  functions  at  r  =  0,  only  the  m  =  1  mode  can  deflect  fluid  off  the  axis.  Since 
this  azimuthal  mode  can  also  stimulate  all  other  azimuthal  modes,  it  provides  an  effective  means  for 
introducing  asymmetry  into  the  flow. 

3.1  Asymmetric  disturbance  growth 

The  minimum  axial  velocity  on  the  axis  is  depicted  in  Fig.  8  for  both  the  axisymmetric  and  asymmetric 
solutions  as  a  function  of  inlet  swirl.  These  results  indicate  that  axisymmetric  breakdown  occurs  at 
Q;,  =  1.48+.  The  small  difference  between  this  result  and  the  previous  axisymmetric  Trial  1  is  due  to 
viscous  effects,  as  confirmed  by  repeating  Trial  1  at  the  lower  Reynolds  number. 

The  degree  that  the  perturbation  affects  the  flow  can  be  evaluated  by  considering  the  integrated 
mode  magnitude, 


(6) 
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Figure  8:  Variation  of  minimum  axial  velocity  at  r  =  0  and  Qm  versus  inlet  swirl.  For  asymmetric 
flows,  a  1%  m  =  1  axial  velocity  perturbation  is  introduced. 


where  V  is  a  radial  slice  through  the  nozzle.  For  the  flows  with  axisymmetric  inlet  swirls  less  than  that 
leading  to  breakdown,  an  asymmetric  perturbation  of  the  form  given  by  Eq.  5  does  not  show  significant 
amplification  and  decays  with  axial  distance.  This  behavior  is  reflected  in  Fig.  8  which  depicts  the 
variation  in  Qm  with  inlet  swirl  for  the  converged  asymmetric  solution.  For  a  non-swirling  flow,  the 
disturbance  propagates  directly  down  the  pipe.  The  addition  of  swirl  both  causes  the  disturbance 
to  spiral  as  it  travels  downstream  and  shears  the  perturbation,  increasing  the  dissipation  rate.  For 
these  solutions,  the  weak  asymmetry  gradually  decays  and  does  not  have  any  significant  effect  on  the 
axisymmetric  flow.  As  demonstrated  by  the  small  values  of  Qrn  and  by  the  agreement  between  the 
axisymmetric  and  asymmetric  solution  curves  in  Fig.  8,  significant  asymmetries  do  not  develop  until 
an  inlet  swirl  of  1.49. 

At  =  1.49,  the  axisymmetric  flow  has  just  undergone  vortex  breakdown  and  a  stagnation  point 
has  developed  on  the  axis.  In  the  cases  predicted  to  have  axisymmetric  breakdown,  the  perturbation 
is  amplified  by  several  orders  of  magnitude,  as  depicted  in  Fig.  8,  and  leads  to  large  changes  in  the 
flow  field.  These  changes  are  also  immediately  evident  in  minimum  axial  velocity  which  indicates 
that  the  steady  asymmetric  solution  does  not  develop  a  stagnation  point  on  the  axis  even  after  the 
axisymmetric  solution  is  predicting  breakdown.  While  asymmetric  solutions  have  been  calculated  in 
this  geometry  for  swirl  ratios  up  to  1.8,  these  higher  swirl  cases  do  not  converge  to  a  steady  solution, 
but  continue  to  have  large,  nearly  periodic,  oscillations  for  long  non-dimensional  times. 

Evolution  equations  can  be  derived  for  the  integrated  mode  magnitude  Qrn  by  integrating  the 
Navier-Stokes  equations  over  the  radial  slice  V. 
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where  qm  =  uTO  •  uTO,  Vmf(r,z)  =  eim0V  (f(r,  z)e~tm0) ,  and  :  is  the  dyadic  double  dot  product  op¬ 
erator.  The  left  hand  side  of  each  equation  represents  the  sum  of  the  growth  of  the  mode  magnitude 
within  the  control  volume  and  boundary  effects.  The  right  hand  side  contains  the  source  terms  re¬ 
sulting  from  viscous  dissipation,  non-linear  interactions  among  perturbation  modes,  and  focusing  of 
perturbations  by  the  base  strain  rate[37].  For  each  mode,  the  viscous  dissipation  always  leads  to  a 
decrease  in  mode  magnitude.  The  non-linear  interactions  of  each  mode  can  be  identically  matched 
to  terms  with  opposite  sign  in  other  modes.  Thus,  while  energy  is  transferred  between  modes,  this 
mechanism  does  not  lead  to  an  overall  increase  in  the  strength  of  the  perturbations.  For  the  axisym- 
metric  mode,  this  term  can  be  related  to  the  work  required  to  maintain  the  asymmetry  and  tends 
to  initially  decrease  the  axisymmetric  velocities.  This  energy  transfer  corresponds  to  the  additional 
“dissipation”  used  by  Leibovich  and  Kribus[16]  to  stabilize  their  axisymmetric  soliton  solution  at  large 
enough  Reynolds  numbers  that  viscosity  alone  was  insufficient. 

The  remaining  volume  source  term, 

i  i  /  u„ ,  um  .  V0U0  dV, 

Jv 

is  the  only  term  for  each  Fourier  mode  that  can  change  the  total  magnitude  represented  by  the  sum  of 
all  Qrn .  This  term  represents  a  focusing  of  asymmetric  modes  by  the  axisymmetric  strain  rate  of  the 
base  flow  and  leads  to  perturbation  growth  in  regions  where  the  base  flow  tends  to  be  compressing  in 
the  direction  of  the  perturbation  velocity.  Thus,  for  near-columnar  base  flows  which  possess  low  strain 
rates,  the  perturbation  growth  is  small.  However,  in  flows  containing  breakdown,  the  region  near  the 
leading  edge  of  the  recirculation  zone  has  large  strain  rates  and  decelerated  flow  both  of  which  tend 
to  amplify  asymmetric  disturbances. 


3.2  Asymmetric  flow  evolution 

In  this  section,  we  examine  the  transient  evolution  of  flow  asymmetries  on  initially  axisymmetric, 
swirling  base  flows.  In  particular,  we  show  that  for  asymmetric  modes  can  provide  a  relieving  effect 
such  that  non-breakdown  solutions  can  be  extended  to  slightly  higher  swirl  ratios  (i.e.  H  >  fij). 
Furthermore,  we  show  how  bubble  modes  of  breakdown  transform  into  spiral  modes  with  the  presence 
of  asymmetries. 

To  examine  the  onset  of  breakdown  in  a  flow  with  asymmetries,  simulations  are  performed  in 
Pipe  A  with  inlet  swirls  larger  than  f Figures  9  and  10  depict  the  evolution  of  the  flow  after  the 
perturbation  is  introduced  at  the  inlet.  The  axisymmetric  base  flow  has  an  inlet  swirl  of  1.55,  greater 
than  lb,  =  1.48,  and  an  axisymmetric  recirculation  zone  exists  at  t  =  0  in  the  test  section.  Figure  9 
suggests  that  there  are  large  oscillations  in  the  nature  of  the  solution  beginning  at  about  t  =  40  non- 
dimensional  times  {^Woo/S).  At  this  point,  the  asymmetries  have  convected  to  the  breakdown  bubble 
and  are  beginning  to  be  amplified.  This  is  also  reflected  in  the  temporal  evolution  of  the  minimum 
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b.  Minimum  axial  velocity  on  axis. 

Figure  9:  Development  of  integrated  mode  magnitudes  and  minimum  axial  velocity  at  r  =  0  as  a 
function  of  time  for  swirl  ratio  =  1.55. 


axial  velocity  on  the  axis.  This  quantity  shows  large  fluctuations  both  positive  and  negative  as  the 
stagnation  point  moves  off  the  axis,  vanishes,  and  reappears.  These  changes  are  also  clearly  evident  in 
the  instantaneous  streamline  images  of  Fig.  10.  Initially,  the  bubble  is  seen  to  grow,  both  in  velocity 
magnitude  and  size,  with  the  significant  asymmetric  velocity  confined  to  the  recirculation  zone.  The 
asymmetries  continue  to  grow  in  time  such  that  the  stagnation  point  is  deflecting  off  of  the  r  =  0 
axis  leading  to  the  characteristic  corkscrew  of  a  spiral  breakdown  (t  =  60).  For  t  =  80,  the  flow  is 
deflected  less  strongly  and  the  strain  rates  decrease.  This  leads  to  a  decrease  in  the  asymmetries  and 
the  stagnation  point  reappears  near  the  axis,  leading  again  to  increases  in  the  asymmetric  nature  of 
the  flow.  As  the  stagnation  point  returns  to  the  axis,  the  helix  angle  of  the  spiral  is  increased  and 
the  nearly  axisymmetric  bubble  mode  reappears  (t  =  110).  This  cycle  continues  with  diminishing 
amplitude  for  a  considerable  time.  Eventually,  the  solution  converges  to  an  approximately  steady  flow 
with  no  stagnation  point.  While  asymmetry  remains,  it  does  not  amplify  significantly  because  of  the 
much  smaller  strain  rate  in  the  axisymmetric  component  of  the  solution. 

This  three-dimensional  solution  closely  resembles  a  continuation  of  the  columnar  solution  branch 
that  exists  for  purely  axisymmetric  flows  just  before  breakdown.  Thus,  while  Wang  and  Rusak[38] 
have  demonstrated  that  axisymmetric  flows  with  swirl  ratios  greater  than  the  critical  swirl  ratio 
(i.e.  subcritical  flows)  are  unstable  to  axisymmetric  disturbances,  this  result  indicates  that  these 
flows  may  be  stabilized  by  the  introduction  of  asymmetric  components.  This  result  has  also  been 
obtained  by  solutions  in  different  pipe  geometries  with  a  variety  of  perturbation  magnitudes  less  than 
5  percent  of  the  base  flow[37].  It  appears  that  the  additional  degrees  of  freedom  created  by  removing 
the  axisymmetric  constraint  provides  a  relieving  effect  that  allows  the  non-breakdown  solution  to  be 
extended  to  at  least  slightly  higher  swirl  ratios.  This  is  consistent  with  the  consideration  of  weak 
asymmetries  as  a  “dissipation” ,  provided  that  the  mean  strain  rates  are  small  enough  that  significant 
amplification  of  the  perturbations  does  not  occur. 

Simulations  at  a  higher  swirl  ratio  of  1.70  show  a  behavior  similar  to  the  evolution  described  above, 
but  after  400  non-dimensional  time  units,  a  periodic  solution  exists  with  no  indication  of  damping. 
This  solution  oscillates  between  a  slightly  asymmetric  recirculation  bubble  and  a  cork-screw  vortex 
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b.  t  =  60 


f.  t  =  170 

Figure  10:  Instantaneous  streamlines  of  flow  after  asymmetric  perturbation  at  an  inlet  swirl  of  1.55 
in  Pipe  A  with  Re  =  500. 


core  following  the  mechanics  described  above.  Seven  cycles  with  a  period  of  about  37  non-dimensional 
time  units  were  simulated  with  no  appreciable  convergence  to  a  specific  solution. 

4  Summary 

We  have  found  that  the  onset  of  axisymmetric  vortex  breakdown  in  swirling  pipe  flows  occurs  when 
the  vortex  attains  local  criticality.  Criticality  was  determined  by  solving  a  local  eigenvalue  problem 
for  the  wavenumber  of  standing  waves  on  the  local  mean  flow.  Approximate  determination  of  flow 
criticality  using  the  local  swirl  number,  while  providing  correct  qualitative  information,  did  not  result 
in  a  unique  maximum  swirl  number  above  which  vortex  breakdown  occurs.  A  transient  simulation  of 
the  evolution  of  vortex  breakdown  revealed  that  downstream-running  waves  are  trapped  approximately 
at  the  location  of  flow  criticality.  These  trapped  waves  are  slowly  amplified  and  eventually  result  in  a 
breakdown  region  with  reversed  flow  and  enlarged  core  size. 

Non-axisymmetric  disturbance  were  found  to  decay  on  base  columnar  base  flows;  however,  for  base 
flows  with  bubble  breakdown,  non-axisymmetric  disturbances  were  amplified  as  a  consequence  of  the 
larger  strain  rates  that  were  present  in  the  initial  flow.  The  amplified  disturbances  resulted  in  the 
formation  of  spiral  breakdown  if  the  inlet  swirl  were  sufficiently  large.  For  a  slightly  smaller  inlet  swirl 
that  would  be  predicted  to  undergo  axisymmetric  breakdown,  near  columnar  solutions  were  found. 
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This  suggests  that  the  energy  transfer  to  the  asymmetric  modes  acts  as  a  stabilizing  influence  until  a 
sufficiently  large  amplitude  is  attained. 
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Author:  Dr.  Darmofal 

Question  by  Dr.  Greenwell:  What  is  the  effect  of  the  tube  wall  constraint  on  the 
structure  (and  stability)  of  the  breakdown.  Specifically,  do  stable,  closed  bubbles  form  in 
unconstrained  vortex  flows? 

Answer:  At  this  point,  we  do  not  know.  For  the  pipe  flow,  we  have  found  that  spiral 
breakdown  is  the  final  form  for  all  of  the  simulations  we  have  performed.  But,  even  this 
conclusion  is  based  on  swirls  that  are  only  slightly  (10-20%)  higher  than  the  critical  swirl 
for  axisymmetric  breakdown. 

Question  by  Dr.  Mitchell:  Have  you  studied  other  vortex  models  beyond  Burgers’ 
vortex  model?  With  delta  wings,  experiments  have  shown  a  trend  to  have  a  vortex 
different  from  that  demonstrated  by  Burgers’  vortex. 

Answer:  Yes.  We  have  looked  a  less  physical  models  (such  as  the  Rankine  vortex  with  a 
plug  flow  axial  velocity)  as  well  as  profiles  from  experimental  data.  The  conclusions  are 
quite  similar. 
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